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Summary. The social- force model is systematically modified to achieve a satisfy- 
ing agreement with the fundamental diagram. Furthermore, our modification allows 
an efficient computation of the simulation. Finally, different simulation-results will 
be compared to empirical data. 



1 Introduction 

Microscopic models are state of the art for computer simulation of pedestrian 
dynamics. The modeling of the individual movement of pedestrians is used 
to design a description of macroscopic pedestrian flow which allows e.g. the 
evaluation of escape routes, the design of pedestrian facilities and the study 
of more theoretical questions. For a first overview see [1-3]. 

One primary test, whether a model is appropriate for a quantitative de- 
scription of pedestrian flow, is the comparison with the fundamental diagram 
[4-7]. 



2 Social-force models 
2.1 General considerations 

The social- force model was introduced in [8] . It is expressed by the equation 
of motion 

TOi ^i = m ~ = F t = Ff rv + F[ ep + fluctuations. (1) 

where x$ is the position of the i th pedestrian, v\ its velocity and m, its 
mass. 

n 

^r p = E^ ep (w^),^N (2) 
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is the force due to interaction with other pedestrians. In the original model 
the introduction of the repulsive force F[ ep between pedestrians is motivated 
by the observation that pedestrians try to keep distance to others to secure 
their 'private sphere'. This behavior is also observed for the environment, i.e. 
pedestrians do not walk too close to walls and stairs. Ff rv is a driving term 
that makes the pedestrians attempt to move with their intended velocity. 

The term fluctuations is added to take into account random variations of 
the behavior of pedestrians. Additionally fluctuations arise from deviations 
from the usual rules of motion. This term becomes important in high density 
situations [8]. 



2.2 Simplified one-dimentional realization 

For the modeling of pedestrian dynamics the pedestrians are represented by 
one-dimensional 'circles' with velocity dependent length di and position (of 
the center of mass) Xj, moving in a continuous space. 

An important aspect is the dependency between the current velocity and 
the required space of pedestrians [9]. As mentioned in [10] it is possible, on 
basis of empirical results, to determine the required length di for the i th 
pedestrian, that moves in a single file with a velocity m, with 0.1 m/s < Vi < 
l.Om/s: 

di = a + bvi with a = 0.36 m and b— 1.06 s (3) 

One fundamental quantity in the model is the distance between the i th 
pedestrian and the j th pedestrian which is defined as: 

dist i , j (t)=A(x i , j )-^±^- (4) 

with A(xij) — \xi — Xj\ the distance between the center of mass of the i th 
pedestrian and the j th pedestrian. |.| denotes the absolute value in K. 




L, — dial — J 

Fig. 1. Illustration of two pedestrians 



The velocity Vi and the centre of mass Xi are calculated by integrating 
Eq. (1) once and twice, respectively. 
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For simplicity fluctuations in (1) are neglected and thus the forces are 
reduced to a driving and a repulsive term = Ff rv + F^ ep . It is assumed 
that in the one-dimensional case a pedestrian only interacts with the man 
directly in front. Consequently (2) reduces to 

Fr = F reP i{x . +uXhVi) (5) 

According to [8] we choose the driving term 

F frv =m .^i i (6) 
T 

where is the intended velocity of the i th pedestrian and r is the ac- 
celeration time. For simplicity we set > 0, x i+ i > Xi and to* = 1. The 
repulsive force is modeled as: 

*cond.: {disti, i+1 (t) < aA j( ";- H> < 0)V (dist iti+1 (t) > / j A d(d "^-+ l) > 0). 

The quantities a and /3 are chosen such that a < (3. To avoide the case 
where two pedestrians pass each other, we assume that the i th pedestrian 
stops once she/he is in contact with the pedestrian in front, i.e. disti t i+i = 0. 



3 Motivation for Event-Driven Simulation 

An event is a single occurrence of a change in the forces acting on a pedestrian 
e.g. when the i th pedestrian comes 'to close' to the pedestrian in front, the 
force Fi changes at a special time - t event - from Ff rv to Ff" v + F[ ep . 
The search for eventually occurring events is performed stepwise with a time 
bound At. It should not be so small as to make the simulation run too long, 
nor should it be so large as to make the number of events unmanageable. Since 
the forces acting on the pedestrians do not change between two events (eventl 
and event2), on can integrate Eq. (1) from t event \ to t even t2 to calculate the 
velocity and the position of the pedestrians in the interval [t even ti t even t2\- 

This approach uses a list of events that occur at various times, and handles 
them in order of increasing time. The simulation makes time 'jump' to the 
time of the next event. The simulation proceeds event-by-event rather than 
step-by-step. 



4 Event-driven simulation with velocity-adaptation 



The pedestrians shall adjust their actual velocity in such a way that they 
keep a minimum distance to the pedestrian in front. Again the forces are 
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only influenced by actions in front of a pedestrian. At the moment where 
the distance between the i th pedestrian and the pedestrian in front reaches 
decreasingly a, the i th pedestrian is subject to the repulsive force (7). This 
leeds to a reduction in the velocity u,. Once the distance between the i th 
pedestrian and the pedestrian in front reaches increasingly (3 the action of 
(7) vanishes. The velocity of the i th pedestrian increases again. 



velocity2 




t_decelerafion t_acce I e ration 



Fig. 2. The variation of the velocity of a pedestrian 

As the required length d is modeled as a linear function of the velocity, the 
quantities a and (3 are chosen to be linearly velocity-depended. (3 is defined 
as the difference between the required length d% and the step-length stj(t) of 
the i th pedestrian and a as See Fig. (4). 

Because of 

sti = st a + stbVi with st a — 0.235 m and s£& = 0.302 s (8) 

according to Ref. [11] and Eq. (3) one gets the complete definition of 
= (3 a + v {(3 a = 0.125 m, (3b = 0.758 s) and thus of a = a a + abv 
(a a = 0.0625m, a b = 0.392s). 

From the definition of the repulsive force (7) one recognizes that times, 
at which the distance between two pedestrians is , a or (3, are times where 
events occur, i.e. at those times force-variations happen. To calculate those 
events, one solves the equation: 

disU{t) = C,C G {0, a, (3} (9) 

5 Results 

To compare the model result with the empirical fundamental diagram of 
the single-file movement [10], a periodic passageway of length L = 17.3 m is 
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Fig. 3. Illustration of a, and d 



selected. The values for the intended velocity v° is distributed according to 
a normal-distribution with a mean value of /i = 1.24m/s and a — 0.05 m/s . 
According to Helbing [8], t = 0.61s is a reliable value. 

For every run at t = all velocities are set to zero and the pedestrians 
are uniformly distributed. After 300 s relaxation, measurement are performed 
300 s long. The determination of velocity mean value over all pedestrians is 
done after every event-evalutation. 



1.4 



L = 1 7.3m, 0a= 0.1 25m, |3b = 0.758s 
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Fig. 4. Results: The fundamental diagram 



Figure (5) presents the dependency between mean velocity and density 
for different approaches to the velocity function (with and without adapta- 
tion) . To demonstrate the influence of the velocity dependence of the required 
length different values for the parameter b were selected. 

With b = 1.06 s and the repulsive force (7) a good agreement at low 
densities between the velocity-density relation predicted by the model and 
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the empirical fundamental diagram is obtained. However at high densities 
the velocity-density relation predicted by the model is inaccurate. 

6 Discussion and summary 

From the empirical fundamental diagram one determines b — 1.06 s and a = 
0.36 m, see [10]. The parameters for the step-length st are also empirically 
known, see [11]. Thus the only free parameter of the model is a. The concept 
of velocity adaptation and the empirically known parameters lead to a good 

agreement with the fundamental diagram at low densities if one chooses a = 



2 

However at high densities one observe a discrepancy between the obtained 
results and the fundamental diagram. This can be explained by the fact, that 
at high densities the reaction-time and the behavior of pedestrians become 
more random so that the deterministic character of the model don't reflect 
anymore the stochastic behavior of the pedestrians. An improvement of the 
model should include a stochastic in the reaction-time of pedestrians and 
determining an appropriate ratio of a to (3 (here it was set to \). 
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